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Starting from the luminosity-redshift relation recently given up to second order in the Poisson 
gauge, we calculate the effects of the realistic stochastic background of perturbations of the so- 
called concordance model on the combined light-cone and ensemble average of various functions of 
the luminosity distance, and on their variance, as functions of redshift. We apply a gauge-invariant 
light-cone averaging prescription which is free from infrared and ultraviolet divergences, making 
our results robust with respect to changes of the corresponding cutoffs. Our main conclusions, in 
part already anticipated in a recent letter for the case of a perturbation spectrum computed in 
the linear regime, are that such inhomogeneities not only cannot avoid the need for dark energy, 
but also cannot prevent, in principle, the determination of its parameters down to an accuracy of 
order 10"'^ — 10~^, depending on the averaged observable and on the regime considered for the 
power spectrum. However, taking into account the appropriate corrections arising in the non-linear 
regime, we predict an irreducible scatter of the data approaching the 10% level which, for limited 
statistics, will necessarily limit the attainable precision. The predicted dispersion appears to be 
in good agreement with current observational estimates of the distance-modulus variance due to 
Doppler and lensing effects (at low and high redshifts, respectively), and represents a challenge for 
future precision measurements. 

PACS numbers: 98.80-k, 95.36.-|-x, 98.80.Es 



I. INTRODUCTION 



In a recent letter [1] we have presented the main ideas and most significant results of a preliminary study of the 
effects of a stochastic background of inhomogeneities on the determination of the dark-energy parameters in the 
context of modern precision cosmology. The main conclusions of that analysis, based on the use of a perturbation 
spectrum valid in the linear regime, were as follows. On the one hand, such kind of perturbations cannot simulate 
a substantial fraction of dark energy: their contribution to the averaged flux-redshift relation is both too small (at 
large values of the redshift z) and has the wrong z-dependence. On the other hand, stochastic fluctuations add a 
new and relatively important dispersion with respect to the prediction of the homogeneous and isotropic Friedmann- 
Lemaitre-Robertson- Walker (FLRW) cosmology. This dispersion is independent of the experimental apparatus, of 
the observational procedure, of the intrinsic fluctuations in absolute luminosity, and may prevent a determination of 
the dark-energy parameter flAiz) down to the percent level - at least if we are using the luminosity-redshift relation 
alone. Another important conclusion was that (light-cone averages of) different functions of the same observable 
get biased in different ways, with the energy flux sticking out as the observable which gets minimally affected by 
inhomogeneities, irrespectively of the redshift binning utilized^. 

The method we have followed, in order to arrive at the above-mentioned conclusions, consists of two different 
steps. In the first step we start from an exact expression for the luminosity-redshift relation in the special "geodesic 
light-cone" gauge introduced in [3] . We then transform this expression, up to second order in perturbation theory, to 
another gauge in which perturbations are known up to that order, the so-called Poisson gauge (PG) (see e.g. [4]). The 



Also redshift binning reduces biases and selects the flux as the preferred variable [2]. 
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second step consists of performing the relevant light-cone and ensemble averages, and in inserting a realistic power 
spectrum of stochastic perturbations. The light-cone average procedure appropriate to this context was formulated 
and discussed in [3] , by extending to null hypersurfaces the gauge-invariant procedure for space- like domains previously 
defined in [5, 6] and also applied in [7]. 

Details on the first stage of this two-step process have been presented in a recent paper [8], while in this work 
we provide a detailed implementation of the second step. The computation method is basically the same as the one 
already used in [9] , but it will involve the full second-order results obtained in [8] . We will consider in detail both a Cold 
Dark Matter (CDM) model and a ACDM one (the so-called concordance model). We will reproduce, in particular, 
the results reported in [1] based on the use of the power spectrum of [10], valid in the linear perturbative regime. 
However, we will also extend our treatment by adding the effects of baryons and by considering two parametrizations 
of the HaloFit model [11, 12], describing the density power spectrum in the non-linear regime. 

The paper is organized as follows. In Sect. II we recall, for the sake of completeness, the results of [8] for the 
contribution of scalar perturbations to the light-cone average of the flux-redshift relation, to second order, in the 
Poisson gauge. We also reorganize the many different contributions in a convenient form for the actual estimates 
to be carried out. In Sect. Ill we present a few important aspects and consequences of the process of combining 
light-cone and ensemble averages, considering in particular the luminosity distance and its phenomenologically most 
relevant functions. We also introduce a convenient spectral parametrization of the inhomogeneous averaged terms. In 
Sect. IV we discuss some relevant results about the dynamical evolution of scalar perturbations, up to second order, 
required for the computation of their averaged contribution. In Sect. V, as a warm-up exercise, we apply our methods 
to a simple perturbed CDM model, where most calculations - except for the explicit mode integration over the given 
power spectrum - can be done analytically. In the relevant range of z, the full leading result can be written in terms 
of an explicit (and simple) function of z times a particular moment of the spectrum, and this shows that such a model 
badly fails in explaining the data, both in magnitude (in particular, at large redshift) and in z-dependence. 

We then turn to the case of a perturbed ACDM model, where calculations are more involved. They are simplified 
by restricting our attention to the so-called "enhanced terms" (i.e the dominant ones in the relevant range of z), 
already identified in the CDM case. In order to discuss a realistic perturbation background, we will also consider 
a power spectrum which includes the contribution of baryonic matter and takes into account the "Silk-damping" 
effect. In Sect. VI we restrict our computations to the linear power-spectrum proposed in [10], and we evaluate the 
impact of stochastic inhomogeneities not only on the averaged flux-redshift relation, but also on other functions of 
the luminosity distance (in particular, on the distance modulus used in the analyses of the Supernovae data). We 
also discuss the dispersion induced by the presence of the perturbations. In Sect. VII we take into account the effects 
of the non-linear regime by using the HaloFit parametrizations of [11, 12]. We also compute the variance/dispersion 
expected in the distance modulus and attempt a first comparison with available data and phenomenological fits, in 
particular for what concerns lensing at large rcdshifts. Finally, in Sect. VIII we summarize our result and draw some 
conclusions. In addition, in Appendix A we discuss why vector and tensor perturbations, although interesting on 
their own, do not contribute to the averages discussed in this paper. In Appendix B we report the explicit results for 
various spectral coefficients used in the analysis of the CDM model. 

II. AVERAGING THE LUMINOSITY FLUX AT SECOND ORDER 

In this section we recall previous results on the light-cone average of the luminosity flux, ($) ^ {d^) (where 
is the luminosity distance), first computed in a generally inhomogeneous metric background, and then specialized to 
the case of a spatially flat FLRW metric perturbed, to second order, by the presence of small fluctuations of scalar, 
vector and tensor type. 

A. Exact expression for (d^^) in the geodesic light-cone gauge 

When describing the propagation of light emitted by sources lying on the past light-cone of a given observer, it 
is convenient to identify the null hypersurfaces along which the photons reach the observer with those on which a 
null coordinate takes constant values. For this reason we have introduced in [3] an adapted system of coordinates - 
defining what we have called the "geodesic light-cone" (GLC) gauge - in which several quantities greatly simplify, while 
keeping all the required degrees of freedom for applications to general geometries. The coordinates = {t,w,6°') 
(with a = 1,2, 0^ = 9, 0"^ = (j>), specifying the metric in the GLC gauge, correspond to a complete gauge fixing of the 
so-called observational coordinates, defined e.g. in [13, 14] (see also [15]). 

The GLC metric depends indeed on six arbitrary functions (a function T, a two-dimensional "vector" U°' and a 
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symmetric matrix 7afc), and its line-element takes the form^: 

d4^P = T2dw2_2TrfwdT + 7„6(rfr-[/°dw)(d^''-C/''dw) . (2.1) 
In matrix form, the metric and its inverse are then given by: 



-T \ / -1 -^VT 

^GLC^ I _^ T^2^^2 g%\c^\ -T-^ 1 , (2.2) 



where = (0,0), Ub = {Ui,U2), while the 2x2 matrices 7a6 and 7"'' lower and raise the two-dimensional indices. 
Clearly ?z; is a null coordinate (i.e. d^wd'^w = 0), and a past light-cone hypersurface is specified by the condition 
w — const. We also note that ~ defines a geodesic flow, i.e. that (9'^t) Vi, [d^r) — (as a consequence of 
g'^'^ = — 1). Such a 4- velocity defines geodetic observers corresponding to the static ones in the synchronous gauge [9]. 
Let us also remark that, in GLC coordinates, the null geodesies connecting sources and observer are characterized by 
the simple tangent vector /c^ = g'^'^d^w — g'"" — —Sf^T~^, meaning that photons travel at constant w and 9°'. This 
makes the calculation of the redshift and of the area distance particularly easy in this gauge. 

In fact, let us denote by the subscripts o and s, respectively, a quantity evaluated at the observer and source 
space-time position, and consider a light ray emitted by a geodetic source (with four-velocity = — 9^t) lying at 
the intersection between the past light-cone of a given geodetic observer (defined by the equation w — Wo) and the 
spatial hypersurface r = Ts, where for the moment is a constant parameter. The light ray will be received by our 
static geodetic observer at t = Tq > r^.. The redshift Zs associated with this light ray is then given by [3]: 

We shall be interested in averaging the luminosity at fixed redshift, hence on the two-dimensional surface S(wo, z^) 
(topologically a sphere) which lies on our past light-cone {w = Wo) and is associated with a fixed redshift [z = Zg). In 
terms of the r coordinate such a surface corresponds to the equation r = Ts{wo, Zs,6°-) enforcing Eq. (2.3). Hereafter 
Ts will denote this (in general angle-dependent) quantity. 

Also the area distance dA, related to the luminosity distance of a source at redshift Zg by one of Etherington's 
relations [16]: 

dA = il + Zsr^dL , (2.4) 
takes a particularly simple form in the GLC gauge. A direct derivation [3] starts from its general definition [17]: 

< - 1 ■ 

where dflo is the infinitesimal solid angle subtended by the source at the observer position, and dSs is the area element 
on the surface orthogonal to both the photon momentum and to the source 4-velocity at the source's position. It is 
easy to check that the surface dSs is characterized by having constant w and r, and that the induced 2-metric on it 
is nothing but jab [8]. Therefore: 

4^^^^, (2.6) 

d^esine sine 

where we have used the fact that photons travel at constant 0°. Our averaging surface S(wo, Zs), being one of constant 
z, differs from the one of constant w and r but - amusingly - the same formula holds, locally, for the area element 
on it, so that Eq. (2.6) can also be used for our light-cone averages [8]. 

The above result singles out the received luminosity flux, $ ~ dj^^ — (1-1- z^)^^^^^, as an important - and extremely 
simple - observable to average over the 2-sphere S(wo, Zg) embedded in the light-cone. In fact (see [8] for more details): 

{dl'){w^, zg) = (1 + zg)-'l^^ = (1 + zg)-^l^ = (1 + zs)-'^^^, , (2.7) 

J dS J dS A[Wo, Zs) 



^ We have put tildas on the GLC gauge 9° coordinates in order to be consistent with our previous notations in [1, 8, 9]. 
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where 



A{Wo, Zs 



(2.8) 



'S(li>o,Zs) 

is the proper area of S(wo, Zg) computed with the metric ^ab^ and expressed in terms of internal coordinates (wq, z^) 
parametrizing the deformed 2-sphere ^(tUojZs). 

Eq. (2.7) holds non-perturbatively for any space-time geometry, and is the starting point for the computation of 
the average flux summarized in [1] and presented in details here. By using the notations introduced in [1] we can 
write, in particular. 



(fi^^)(Wo,Zs) = {l+ZsY 



ana 



At: 



-^{wo,T,{w,,z,,r),r) 



where we have defined 



Itl,{Wo,Zs) 



A{Wo,Zs) 



47r a(?7i"^)Ar7 



(2.9) 



(2.10) 



and where d^^^^ 



(1 + ZsYa{rls'^)/S.ri is the luminosity distance for the unperturbed FLRW geometry, with scale 
factor aijf). Here rj is the conformal time coordinate, At^ — ijo — Vs^K and we have denoted with rii^^ the background 
solution of the equation for the source's conformal time 77^ = risizsjO"") (see [8, 9]). Note that, according to the above 
equation, the interpretation of I^{wo, Zg) is straightforward: it is simply the ratio of the area of the 2-sphere at redshift 
Zs on the past light-cone (deformed by inhomogeneities) , over the area of the corresponding homogeneous 2-sphere. 



B. Second-order expression for (d^'^) in the Poisson gauge 

Let us consider a space-time geometry that can be approximated by a spatially flat FLRW metric distorted by the 
presence of scalar, vector and tensor perturbations. In the so-called Poisson gauge (PG) [4] (a generalization of the 
standard Newtonian gauge beyond first order), the corresponding metric (in Cartesian coordinates) takes the form: 

ds^C = a^i-q) {-(1 + 2^)drf + 2uj,dTjdx' + [(1 - 2*)% -|- h,j] dx'dx^} . (2.11) 

Here $ and 'J are scalar perturbations, uji is a transverse vector perturbation {d'^uJi = 0) and hij is a transverse and 
traceless tensor perturbation {d^hij = = hi). This metric depends on six arbitrary functions, hence it is completely 
gauge fixed. By including first-order and second-order contributions, the (generalized) Bardeen potentials and ^ 
can be defined as follows: 

$ = ^ + 1^(2) ^ VI/ = ^ + 1^-^'' , (2.12) 

where we have assumed the absence of anisotropic stress in order to set — ^ — ip at first order. 

It is important to stress, at this point, that for the purpose of this paper we can safely restrict our subsequent 
discussion to the case of pure scalar perturbations. In fact, it is true that at second order different perturbations get 
mixed: vector and tensor perturbations are automatically generated from scalar perturbations (see e.g. [18, 19]), while 
second-order scalar perturbations are generated from first-order vector and tensor perturbations. However, a single 
vector or tensor perturbation does not contribute to our angular averages on Y,{'Wo,Zs) (see Appendix A for further 
details). Furthermore, we will treat uJi and hij as second order quantities. In other words, we assume the first-order 
perturbed metric to be dominated by scalar contributions (which is indeed the case if perturbations are generated by 
a phase of standard slow-roll inflation, see e.g. [20, 21]). As a result, we shall also neglect the contributions induced, 
at second order, by first-order vector and tensor perturbations. 

We have already established in [8] a connection between the second-order perturbative expression of the luminosity 
distance and of the integrand of (controlling (d^^)), both written in terms of the PG perturbations and of the 
observer's angles 0°" (remember that photons reach the observer traveling at constant 6*°). In particular, by defining^ 

ntfi - ^i+si^\z.,n+s^i\z.,n , (2.13) 



^ Note the simplified notation with respect to [8]. We have also omitted the indication that w = Wo- 
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we have found that [8] 



and 

M^s)^ |^sm^(l+Ii+Ii,i+l2) , (2.14) 

= 2(5^ + (t. d.)('\ Ii,i+l2 = 24'' + (4^V + (t. d.)('', (2.15) 

where the (t. d.)^^'^^ denote total derivatives terms w.r.t. the angles, giving vanishing contributions either by 
periodicity in (j) or by the vanishing of the integrand at 6* = 0, tt. We only recall here, for later use, the first-order 
total derivative: 

1 rVo _ (0) _ 
(t. d.)W=2J« , 4'^-^/ d7y^^A2^(77,,7„~,7,r), (2.16) 

where A2 = 9~ + cot 6d^ + sin~^ 9d~. It is also convenient to rewrite the PG metric using spherical coordinates (but 
still considering that photons travel at constant 6*"), and define the following quantities [8]: 

Es = l--^ , J={[d+Q]s~[d+Ql)~{[drP]s'[drP]o) . (2.17) 

Here H — dlna/drj, and the lower limit r]in represents an early enough time when the perturbation (or better the 
integrand) was negligible. When ambiguities may occur the superscript (0) denotes the background solution of a 
given quantity (similarly, the superscripts (1), (2) will denote, respectively, the first- and second-order perturbed 
values of that quantity). Finally, we have introduced the useful (zeroth-order) light-cone variables ri± = rj ± r, with 
corresponding partial derivatives: 

d _ 1 



d^^d++d. , dr=d+-d- , d±^ — ^-{d^±dr) . (2.18) 



Also, we shall use hereafter a hat to denote a quantity expressed in terms of the (r/+, 6*°) variables, so that, for 
instance, ii{j{r]^,ri_,6'^) = 'ip{r],r,9°'). 

The results obtained in [8] can then be reported in the following form: 



3 23 



i—l i—1 i — 1 



where Ii, Ii.i, I2 are, respectively, the first-order, quadratic first-order, and genuine second-order contributions of 
our stochastic fluctuations, and where: 

r/^) - -2V(r;W,r("),r) ; T^'^ = 2S,J ; T^'^ = -^Q, ; (2.20) 

r/''^^ = S4V.' - V'o'] ; =Es{{[drP]sr-{[drP]or) i T^'^ = {A + [d+Q]s) [drP]s ; 

ri'''^ = ^S,(7o°^), {2daPsdbPs + daQsdbQs ~ ^daQsdbPs) ; = -5. lim [l^^daPdbP] ; 

Z r— >0 

Tg^'-'^ = 2~sQs (^2aW^ + 2a^Vo - driJ, + 2 dr^'dl^P [rj' , r,„ - 77', ^ ) + [d^^Pl^j ; 

T,^'''^ =2Es:^{[d,i:]s+ns[drP]s) ; ^2Es:^[d'^P]s ; 

its tis 



= -S, / ' drj'^^dr + {drPf+jS'daPdbP] {v' , , ^ i 



r/o = r drj'^dr + {drPf + iS'daPdbP] (77', 0, r ) ; 



rA'''^ = 2S, / da; d+ 



(0) 
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'12 



/ 1 ^ 



'l8 



'21 



'23 



1 



-2(J-2^,) 



/ 

-0- 






A77 


J 'I6 



A77 



J 



Ar] 



1 



'17 — 

(1,1) _ 



■(1,1) _ 



da: 



1 d 



cos I 



,(0) + 



9+Q + -7o"^ daQ d,Q 



A, J ' 



(2.21) 



and 



■(2) _ 



■(2) _ ^ 



a('7o) 



0(2) +^(2)1 (^(0)+,x,r) 



■(2) 



(2) 



'7 



■(2) 



A?7 7^(0) 



0(2) +^(2) 



(ry(")+,x,r) . 



(2.22) 



In the above equations — diag(r 2, r ^ sin ^ 9), and for 7^^^'"'^'' we have used the following identity 



(2.23) 

Let us point out, finally, that the last term in Eq.(2.21) corresponds to a total derivative, and thus to a boundary 
contribution, that superficially looks non vanishing. We believe that this is the result of a naive treatment of the 
angular coordinate transformation, which becomes singular near the poles of the 2-sphere. This contribution, indeed, 
has the same form as the irrelevant one coming from an overall 5*0(3) rotation. 



III. COMBINING THE LIGHT-CONE AND ENSEMBLE AVERAGE OF di AND OF ITS FUNCTIONS 

In the cosmological model we are considering, the deviations from the standard FLRW geometry are sourced by 
a stochastic background of primordial perturbations satisfying -0 = 0, ■0^ 7^ 0, where the bar denotes statistical 
(or ensemble) average. Hence, non-trivial effects on the ensemble average of d^, or of a generic function of it, can 
only originate either from quadratic and higher-order perturbative corrections, or from the spectrum of correlation 

functions such as dL{z,9°-)dL{z' ,9'°-) (see [22]). In this paper, rather than considering the ensemble average of d^, we 
shall deal with that of (d^), where the angular brackets refer to the light-cone average defined in [3] and presented 
in Sect. II (see [23-25] for previous attempts of combining ensemble and space-time averages in the case of spacelike 
hypersurfaces) . 

As already stressed in [9], given the covariant (light-cone) average of a perturbed (inhomogeneous) observable 
5, the average of a generic function of this observable differs, in general, from the function of its average, i.e. 
{F{S)) 7^ F{{S)) (as a consequence of the nonlinearity of the averaging process). Expanding the observable to 
second order as 5 = 5*0 -|- 5*1 -I- 5*2, one finds [1] 

WS)) = F{So) + F\So){Si + S2) + F"(5o)(5p), (3.1) 

where in general (Si) 7^ as a consequence of the "induced backreaction" terms, arising from the coupling between 
the inhomogeneity fluctuations of S and those of the integration measure (see [9]). Eq. (3.1) implies, in our case, that 
different functions of the luminosity distance (or of the flux) may be differently affected by the process of averaging 
out the inhomogeneities, and may require different "subtraction" procedures for an unbiased determination of the 
relevant observable quantities. 
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Let us consider, in particular, the luminosity flux $ ~ d^^ (not to be confused with the Bardeen potential!), 
computed in Sect. II. Performing the stochastic average of Eq. (2.9) (and using Eq. (2.14)) we obtain 



where: 



{dl'){z) - (dr^^)-^(/*(z))-i ^ (d™)-2 [1 + Mz)] , (3.2) 



^(z)ee(Ii)2-(Ii,i+I2) . (3.3) 



We can now apply the general result (3.1) to the flux variable, by setting S* = $ and considering two important 
functions of the flux: F(<i>) = $-1/2 ^ ^^^^ ^nd F($) ~ —2. Slogan <i> + const ^ ^ (the distance modulus). They wiU be 
considered in the following sections, together with the flux. For the luminosity distance we can introduce a fractional 
correction fd, in analogy with Eq. (3.2), such that: 



(di)(z) = d™[l + /,(z)] . (3.4) 
Then, by using the general expression (3.1), we flnd: 

u = -lu + lWUW) , (3.5) 



where, in terms of the quantities defined in (2.13), we have ((<I>i/$o) ) ~ Mi^s )^)' ^^'^ where is defined by Eq. 
(3.3). For the distance modulus we obtain, instead. 



~ m''^^''' = -l-25(logio e) [2U - ((<J'i/$o)')J • (3.6) 

We can also consider, for any given averaged variable (S), the associated dispersion as controlling how broad is 
the distribution of a perturbed observable S around its mean value (S). This dispersion is due to both the geometric 
fluctuations of the averaging surface and to the statistical ensemble fluctuations, and is dcflned, in general, by [9]: 



as^^^[S-{S)) ^ = \/(52)-((5)) . (3.7) 
The dispersion associated with the flux is thus given by: 

fT$ - y(($/$o)')-(W*i^)' - , (3.8) 

while for the distance modulus we find: 



= \l (a*') - ((m>)' = 2.5(logio e)v/ ((<i>i/<i>o)') . (3.9) 

The above results will be applied to the case of a realistic background of cosmological perturbations of infiationary 
origin in the following sections. 

Let us conclude this section by introducing a convenient spectral parametrization to be used for the various terms 
contributing to the fractional corrections of our observables, and to the corresponding dispersions. We start by 
recalling that the simplest way to implement the ensemble average of a stochastic background of scalar perturbations 
ip is to consider its Fourier decomposition in the form: 

V'(ry,f) = J d^ke^^^-^Mv)Eik) , (3-10) 

where - assuming that the fiuctuations are statistically homogeneous and isotropic - i? is a unit random variable 

satisfying E*{k) = E{—k), as well as the ensemble-average conditions E{k) = and E{ki)E{k2) = 5{ki + ^2). As a 
simple illustrative example one has: 



cPk d^fc'— 7F— 7^ r (fn r 



(V'.V'.) = / ^^i-^E{k)E{k') I ^ M4°^)e'''''-' ,o> ^fe'(^r)e' 



(0) 



{Q)\ irk'-x 



(0) 

r=rio—7]s 



(27r)3 ^ ' ' 'J 47r 

'^'^ \Mvi'')\'- r^nik,vi'^ (3.11) 



(27r 
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ii's) {4's 



(27r)3 
<fk 



E{k)E{k') 



Ait 



{^)\ irk-x 



(2^) 



d{cos9) ^^f.^^ 
2 



r=rio-V. 
1 



d(cos 0') ikA-n cos e' 
2 



(2^)3 '^^^'^ ^' V fcAry / 



(3.12) 



where in the second hne of both terms we made use of isotropy (i.e. Vfc only dependent on k ~ |fc|), and defined 9 
and 9' as the angles between k and x = rx and between k' and x' = rx' . We have also introduced the (so-called 
dimensionless) power spectrum of '0: 



To give a slightly more involved example, we can consider the average of a term like ^^Qs/ ^W- 

d^kd^k' 



(3.13) 



At] y (27r) 



3 E{k)E{k' 



d-'n 



-2 



V^{k,rio) 



An 

SinInt(fcA?7) 
A:A?7 ' 



V'fe(ryi°^)e^'-"^'' / dryVfe'(?7)e^''-*'''°"") 



(3.14) 



0^ 

A?7 



2 

A?7 



E{k)E{k') 



d^n 



(27r)3 ''"'[y 47r 

°° dA: , sin(fcA77) SinInt(A:A?7) 

—V^(k,rio 



At: 



dr]ipk'{r])e 



ik' -x' {'Qo — rj) 







kl^-q 



kAf] 



where the arrows refer to the case of a time-independent fluctuation mode and 



f 

Sinlnt(x) = / 
Jo 



-dy 



y 



smy 



(3.15) 



(3.16) 



As one can see, the angular average is making the results completely different in the two cases. We remark that the 
presence of the Sinint function is a direct consequence of the integration over time in Qs ■ Consequently, the non-local 
nature of the backreaction terms is reflected in the form of the corresponding spectral coefRcients. 

Following this approach, all the relevant contributions to the averaged functions of the luminosity redshift relation, 
at second order, can be parameterized in the form: 



dk 



(X) = ^ V^ik,rjo)Cx{k,Vo,vi°^). 
{X'){Y') - / -7'^(fc,ry„)Cx'(A:,r;„,,y(0))Cy,(fc,ry„,,7(0)), 



(3.17) 
(3.18) 



valid for any given model of perturbation spectrum. Here, X' and Y' [X) are first (second) order generic functions 
of {r],r,9°'), and the C are the associated spectral coefficients. In the particularly simple case of a CDM-dominated 
geometry the spectral distribution of sub-horizon scalar perturbations is time-independent {drjipk = 0) and, as we shall 
see later, all the spectral coefficients C can be calculated analytically. We should stress, however, that when performing 
numerical calculations the above integration limits will be replaced by appropriate cut-off values determined by the 
physical range of validity of the considered spectrum. 



IV. DYNAMICAL EVOLUTION OF SCALAR PERTURBATIONS UP TO SECOND-ORDER 



For a full computation of the fractional correction /$ what we need, at this point, is the combined angular and 
ensemble averages of the three basic quantities, Ii, Ii.i and l2- We must evaluate, in particular, the spectral 
coefficients {C_(i)}. {C_(i,i)} and {C_(2)}, related to the terms defined in Eqs. (2.20-2.22) in terms of the first and 
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second-order Bardeen potential. For this purpose we need to know the dynamical evolution of the scalar fluctuations, 
at first order for the computation of Xi, Ii^i and at second order for l2- 

Let us consider, first of all, a general model with cosmological constant plus dust sources. For the evolution of the 
scalar degrees of freedom in the Poisson gauge we will follow the analysis performed, up to second order, in [18]. In 
a general ACDM model the linear scalar perturbation obeys the evolution equation 



Considering only the growing mode solution we can set 



(4.1) 



(4.2) 



where 5(77) is the so-called "growth-suppression factor", or - more precisely - the least decaying mode solution, and 
i>o is the present value of the gravitational potential. This growth factor can be expressed analytically in terms of 
elliptic functions [26] (see also [27]), and it is well approximated by a simple function of the critical-density parameters 
of non-relativistic matter (fim) and cosmological constant (IIa) as follows: 



ntl^ - r^A + (1 + ^m/2){i + ' 



(4.3) 



Here f/00 represents the value of 3(77) at early enough times when the cosmological constant was negligible, and is fixed 
by the condition g{r]o) — 1. 

The second-order potentials obey a similar evolution equation, containing, however, an appropriate source term. 
Their final expression in terms of ipo has been given in [18] and reads: 



10, 



where 



5'^ - 3 



01 



(4.4) 



(4.5) 



(4.6) 



and where we have introduced the functions BA{r]) = "Hq ^ iKVo) + 3ftm{T]o)/2] ^ Ba{v)7 with yl = 1, 2, 3, 4, and with 
the following definitions: 



Biiv)-- 

Bsiv)- 
and with 



dr~jn'if]){m-^rCiv,v), B2iv) = 2 / dfjn'in) 2(^(77) -1)2 -3 + 30^(77) ^(77,77) 



d77 e(^7)-f- ^(77,77), ^4(77): 



2/,^^ a(?7) 



a{r]o 



g{ri)n{fi) - g{fi)^^n{r]) 
aw 



dfjC{r],fj) 



(4.7) 
(4.8) 

(4.9) 



Here 77„i denotes the time when full matter domination starts [18]. Its precise value is irrelevant since the region of 
integration around 77,„ is strongly suppressed. Finally, a^i is the so-called non-gaussianity parameter (see [18]), which 
approaches unity in the standard infiationary scenario. 

For further use let us now evaluate the ensemble (and angular/light-cone) average of the different operators defined 
in Eq.(4.6), when applied to diipodjipo- Considering first the ensemble average of O^^ diipodjipo, and Fourier-expanding 
ipo, we get (see [28]): 



2(fc • q) 



(4.10) 



10 



(from this point, and up to the end of this section, we will neglect all suffixes "o" present in terms inside the integrals). 
By using the Taylor expansion of ■0*- ^ around q = Q we have: 



{k ■ qf 



(4.11) 



where fc = q= \q\^ and where the latter terms have been obtained by using the Hessian matrix — dk^dk^^'^ = 
dki[{k^ /k)dktl'l\. Combining the last two results, writing the integral over k as 2T:k^dk J^^ dcosa{...), where 
k ■ q~ kq cos a, and integrating over cos a, we obtain: 



(M 



r(5(^)(g)e*« *'' / 2TTk^dk 



16fc 



15 



(4.12) 



Note that the integrand's dependence on the angle 9 between x and q only arises from the exponential term exp(i(f-xr), 
which disappears in the presence of 5^^\q). As a consequence, the angular average has no impact on this particular 
term and we get: 



dfc- 



27r2 



— V'fcSfeV'fc + -r^ ^kOki'k 
15 15 



(4.13) 



Let us note also that, quite generally, dki^k ^ i^k/k, and thus the above term is of the same order as (tpoipo)- 

le p 

dh 



By repeating exactly the same procedure for the O^^ di^odjtpo term, we find: 



(2vr) 



^ 5^3) (^gig-fir / 4jr/s2^/j 



20 



^kdki^l 



(4.14) 



where the last two contributions are vanishing because of the g-integration. We are thus left with the following simple 
result (insensitive to the angular average, as before) 



'^Ol^ diipodjipo^ = Ol^ ditpodjijJo = 
Finally, the last term O^^ ditpodjtpo — iVi/joP is trivial, and gives 



dk 
~k 



-V4,{k,rio) 



(4.15) 



(4.16) 



V. AN ILLUSTRATIVE EXAMPLE: THE CDM MODEL 



In the previous sections we have computed, up to the second perturbative order, the general form of the corrections 
induced by the stochastic fluctuations of the geometry on the averaged luminosity flux, for a spatially flat FLRW 
metric and for a generic spectrum of metric perturbations. We have found that such corrections are controlled by the 
combined angular and ensemble averages of three basic quantities, Ii, and l2- In this section we will evaluate 
such averages for the particular case of the standard CDM model, as a simple illustrative example. The results we 
will obtain, will give us useful information about the dominant terms to be selected also for the discussion of the 
phenomenologically more relevant case, the ACDM model, presented in the forthcoming sections. 



A. The quadratic first-order contributions (Xi)^ and 

Let us now explicitly calculate the spectral coefficients for the first two backreaction terms (induced by the averaging) 
contributing to Eq. (3.3). The genuine second-order term (I2) will be discussed in the next subsection. We start 
considering the first term, (Ii)^. From Eq.(2.20) we obtain: 

3 3 p 00 77 3 3 

(^ = EE('^^'^)('^-^'^)- / y n(fc:'?o) 5]^C^a,(A:,^„,,7.)C^a)(fc,'7o,r/.) (5.1) 

4=1 j = l i=l j = l ' ' 
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TABLE I: The spectral coefficients C^n) for the (t^'^') (T^^^) terms. 



' i 


■i 




() sin feArj 


'2 


-2(l-«^) (l-^)+2(l-^)i^(cosfcAr,-^) 


T-{1) 
'3 


^SmInt(fcA,7) 



(notice that, from now on, the background solution -qf^ will be simply denoted by rjs). For a dust-dominated phase 
the spectral distribution of sub-horizon scalar perturbations is time independent S^V'fc = 0, and the scale factor a{r]) 
can be written as 0(77) — a{r]o)i'i]/Vo)'^- We can also define 



fo,: 



dv—, ; 



(5.2) 



(recall that 7?i„ satisfies, by definition, rjin <^ Vo,s)- All the spectral coefficients of Eq. (5.1) can then be easily 
calculated, and the result is reported in Table I In a similar way, the second contribution to Eq. (3.3) can be 
expressed as: 



■^0 1=1 



(5.3) 



and the explicit form of the spectral coefficients C. 



for the CDM case, is presented in Appendix B. 



Considering all contributions generated by Xi and Ii.i, we find that the dominant contributions arc all contained in 
(Xi^i), and are characterized by spectral coefficients proportional to k'^ (such dominant terms have been emphasized, 
in the Appendix B, by enclosing them in a rectangular box). They correspond, in particular, to the terms 7^^^'^' with 
{i = 2, 4, 5, 7, 8, 12, 14, 18, 20}. Including only such dominant contributions we find that, to leading order. 



■ 23 



n- 



^ ^ Z^±io fc2 



2(^s 3) /s 2 



Lead 



S '— 



f fsk' 



Js 1.2 Jo ,2 

"^3 "^3 



4S, 
3 



where we have defined 



1 - 



(5.4) 



(5.5) 



and we have used the relation "Hs ~ 2/(3/s). Also, we have included into the function fi,i{z) all the z-dependence of 
these leading contributions. After some simple algebra we find: 



10 - 12^/^TI -f 5z (2 -I- \/rTi) 

27(i + z) (-i + yrr^)^ 



(5.6) 



This function (and thus the corresponding backreaction) fiips sign around z* = 0.205 (as illustrated also in Fig. 1). 

It should be noted, finally, that some of the genuine second-order terms contained into I2 are also associated to 
spectral coefficients proportional to k'^ . However, as we shall see in the next subsection, the corresponding contributions 
to Eq. (3.3) turn out to be roughly an order of magnitude smaller than the above ones, because of approximate 
cancellations. 



We take the opportunity to point out two misprints appearing in Table 2 of [9] where the spectral coefficients of both Ai and A2 should 
have the opposite sign. 
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B. The genuine second-order contribution {I2) 



In order to complete the calculation of /$(z) we still have to consider the genuine second-order backrcaction term 
(Z2). In particular, we must evaluate the spectral coefficients {C^(2)}, corresponding to the various contributions 

defined in Eq. (2.22), in terms of the first-order Bardeen potential. By using the results of Sect. IV, and starting from 
Eqs. (4.4) and (4.5), we can easily see that the only possible fc^-enhanced contributions arising from the coefficients 
{C^{2)} (which only contain functions of (p^^^ and -0^^^) should correspond to the term B3{ri)V^'^did^ (d^tpodjipo) + 

Bi{Ti)d^4'odi'ipo- Therefore, we should obtain 0^^) ~ at leading order. 

Let us discuss and estimate all possible contributions for the CDM model we are considering in this section. In this 
simple case l{ri) = 1 and Bi{rf) — B2{rf) — (see Eq.(4.7)). Furthermore, restricting our attention to the standard 
inflationary scenario, we can set a„i — 1. Eqs. (4.4) and (4.5) thus reduce to: 



(5.7) 
(5.8) 



(5.9) 



Finally, we can use (as before) the scale factor 0(77) — a[rio) (il/Vo) , and obtain (according to Eq. (4.8)): 

, , 20 1 , , 2 1 

i?3(^) = ^^ ' ^4(^)^-7^^- 

We are now in the position of evaluating the genuine second-order terms, by exploiting the results given in Eqs. (4. 13), 

(2) (2) 

(4.15) and (4.16). Following the classification of Eq.(2.22) we can see that the first two terms 

c(rro and Cirri 

exactly cancel for the CDM case (while they cancel only at leading order for the case of a ACDM model). For CDM 
we have, in particular. 



c(r/^)) = -c(r,(^)) = -^(ry,2 



(5.10) 



Another interesting simplification concerns the terms 7^*'^'', Tj"^^ and 7^^^\ namely those terms for which the in- 
tegrand contains cJ^V'*'^'' or drcf'^^'^ ■ From our previous results, in particular from Eqs. (4.13), (4.15) and (4.16), it 
is easy to see that the V'*-^'' and contributions are unchanged when one averages over the 2-sphere (i.e. over 
9 by isotropy), since the 6'-dependence is removed by the presence of S'^^^q). On the other hand, the presence of 
the r-derivative brings a further factor \q\ cos 9, and the g-integration gives zero, even before performing the angular 
average. It follows that, for a general model. 



(2)n 







c{t: 



(2)^ 







(5.11) 



(2) 

The contribution of 7^ , on the contrary, is nonvanishing because of the presence of the partial derivative 9^ , acting 
on the BA{r}) coefficients. For the CDM model we find, in particular. 



126 



(5.12) 



Finally, for the last two terms 7^ and T-^^^' we obtain 



-(2) 



(2)^ 



To 



189A77 

The sum of all contributions then leads to: 
7 



126 



45 l^fcP 
16fc ^kdk^l 



8k' Adli^l 
45 l^feP 

45 l^fcP 



^C(7:(^') = -2 



i=l 



11 
2 



f s f o 

A77 



16fc ^pkdk^^ 
15 \^k\ 



15 \i>k\' 



(5.13) 



(5.14) 



All the above spectral coefficients are now to be numerically evaluated by using the power spectrum of the CDM 
model. We can easily check, however, that the leading /c^-contributions of these coefficients are given by: 



7 

E 

.i=l 



- Lead 



^ f2 f3 _ f3 

I f2 _ f2\ _ J_s J s Jo 

2 ^•'^ ■'"'2 At] 



(5.15) 



13 



where: 

~ 1 2-2v/r+i + z(9~27r+^) 

^^(^) = -i89 — (i + .)(ym-i) — • ^'-''^ 

Such second-order contributions turn out to be about one order of magnitude smaller than the leading contributions 
of the squared first-order terms of Sect. VA (as can be easily checked, for instance, by comparing the plots of 
and /2). 



C. Full numerical results for the CDM model 



At this point, in order to perform the numerical computations, we need to insert the explicit form of the power 
spectrum. Limiting ourselves to sub-horizon perturbations we can simply obtain 4'k , for the CDM model, by applying 
an appropriate, time- independent transfer function to the primordial (inflationary) spectral distribution (see e.g. [20]). 

The power spectrum of the Bardeen potential is then given by: 

r4k)^(^^y A'^T^ik) , Al = A(^^y^ ' , (5.17) 

where T{k) is the transfer function which takes into account the sub-horizon evolution of the modes re-entering 
during the radiation-dominated era, and is the primordial power spectrum of curvature perturbations outside the 
horizon. The typical parameters of such a spectrum, namely the amplitude A, the spectral index Ug and the scale k^, 
are determined by the results of recent WMAP observations [29]. In our computations we will use, in particular, the 
following approximate values: 

A = 2.41xl0"\ n^ = 0.96, /cq = 0.002 Mpc"V (5.18) 

Finally, since our main purpose here is to present an illustrative example, it will be enough for our needs to approximate 
T(k) by the effective shape of the transfer function for density perturbations without baryons, namely T{k) — To(fc), 
where [10]: 

^°('^)= Lo-H^°Co(,) ' M'^)-M2e + 1.8,), ^^o(.) = 14.2 + , (5-19) 

and where fccq — 0.07 flmoh^Mpc"^ is the scale corresponding to matter-radiation equality, with h = 
ffo/(100kms"^Mpc~^). 

We can easily check that the above transfer function goes to 1 for fc <C fccq, while it falls like fc~^ log k for k 3> fceq. 
For the numerical estimates we will use h — 0.7 and we will set ao = 1, i^m — 1- In that case we obtain fceq ~ 
0.036 Mpc~^ (see [10]), and we can more precisely define the asymptotic regimes of our transfer function as Tq ~ 1 
for k < lO^^Mpc^i, and Tq - k-^\ogk for k > 2.5 Mpc"^ 

Following the results of the previous subsections we can now set 

r°° dk r 1 

h{z)^ — 7'^(fc,z = 0)[/i,i(fc,z)-h/2(fc,z)J, (5.20) 

where /i^i and /2 are complicated -but known for the CDM case- analytic functions of their arguments. However, as 
already stressed, the leading contributions in the range of z relevant to dark-energy phenomenology are sourced by 
terms of the type /(fc, z) ^ (k/'Ho)'^ fiz). In that range of z we can thus write, to a very good accuracy: 

dk f k ^ ^ 



/$(z)~ ^/ia(z)-f /2(^)jy^ T (,?rj ^^-^^^ 

where fi,i{z) and f2{z) are given, respectively, by Eqs. (5.6) and (5.16). 

To proceed, we need to insert a power spectrum as well as infrared (IR) and ultraviolet (UV) cutoffs in (5.20). The 
former can be identified with the present horizon however, considering the used spectra, larger scales give a 

completely negligible contribution. On the other hand, in spite of the fact that our expressions converge in the UV 
for any reasonable power spectrum, some mild sensitivity to the actual UV cutoff will be shown to occur in certain 
observables. The absolute value (and sign) of f<s>{z) for the CDM model, obtained from both Eq. (5.20) and (5.21), 
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FIG. 1: The fractional correction of Eq. (5.20) (solid curve), compared with the same quantity given to leading order by 
Eq. (5.21) (dashed curve), in the context of an inhomogeneous CDM model. We have used for the spectrum the one defined 
in Eq. (5.17). The plotted curves refer, as an illustrative example, to an UV cutoff kuv = lMpc~^. 

are illustrated in Fig. 1, where we can explicitly check the accuracy of the leading order terms (5.21). The figure also 
confirms that the backreaction of a realistic spectrum of stochastic perturbations induces negligible corrections to the 
averaged flux at large z (the larger corrections at small z, due to "Doppler terms", have been already discussed also 
in [9]). In addition, it shows that such corrections have the wrong z-dependence (in particular, they change sign at 
some z) for simulating even a tiny dark-energy component. 

VI. THE ACDM MODEL: POWER SPECTRUM IN THE LINEAR REGIME 

We will now extend the procedure of the previous section to the case of the so-called "concordance" cosmological 
model, using first a power spectrum computed in the linear regime, and then adding the effects of non-linearities 
following the parametrizations proposed in [11] and [12]. In both cases, we will restrict our attention only to the 
/c^-enhanced terms already identified in the CDM model for the light-cone average of the flux variable, as well as to 
the fc'^-enhanced terms which, as we will see, will appear in the variance, or in the averages of other functions of dL^z). 

Hereafter in all numerical computations we will use, in particular, the following numerical values: Oao — 0.73, 
Qrno = 0.27, Qbo = 0.046, and h = 0.7. 



A. Second-order corrections to the averaged luminosity flux 



The power spectrum of the ACDM model is, in general, time-dependent. Considering for the moment only the 
linear regime, the scalar power spectrum can be written, starting from Eq.(4.2), as 



.givo) 



(6.1) 



and in this case we can easily extend the results previously obtained for the CDM model, concerning the leading 
(/c^-enhanced) contributions to the averaged-flux integral J^. Using the general definitions of the T^^^'^"* terms, we 
first notice that such enhanced contributions arise from the following particular (sub)-terms appearing in Eq. (2.21): 



'2,L 



T 



(1,1) 



5,L 



-S, lim [-f^'daPdbP] 



15 



'8,L 



= 2S 



1 



'12, L 



^2 



K 



{{[drP]sf + {[drP]of 



T 



(1,1) 



14, L 



(1,1) 



18, L 



■{lf')sdaQsdb[drP]. 



T-(l,l) _ oJ_ 
'20, L ^ 



[drP]s[dri^]s 



■2) 



where we have added the suffix "L" to stress that we are reporting here only the leading terms. 

In order to calculate the corresponding spectral coefficients, we note that their time dependence can be factorized 
with respect to the fc-dependence whenever the time variable does not appear in the exponential factor exp(i/c • x) 
present in our integrals (see for instance Eqs. (3.11), (3.12), (3.14), (3.15)). This is indeed the case for the terms 
T2^i^\ tI^I^\ '^^L^^ Til'i^, Till! 3'iid 72q'l ■ In that case the previous CDM results for the leading spectral 

coefficients can be simply generalized to the ACDM case through the following procedure: (i) by inserting a factor 
diVs)/ giVo) whenever ips is present in the initial term; and (ii) by replacing the fo^s factors (see Eq.(5.2)), arising 
from the presence of Pq s terms, by: 



/o,. 



drj 



a'{Vo,s) g{r]o)' 



(6.3) 



For the remaining two terms T^^j^'^ and Ti^'^' the integrals are performed along the path r = ijo — 1], and the 
time dependence cannot be fully factorized. Therefore, the evaluation of the double integrals over ry and k is much 
more involved than in the previous cases. However, a good approximation of the exact result can be obtained by 
replacing 'ip{r]',rjo — rj',6'^), appearing in the integrands of T^^£^^ and Ti^'l\ with ip{r]s,rio — r/'i^"^) (this is so since 
the leading contributions to the time integral arise from a range of values of rj approaching rjg). By adopting such an 
approximation we can follow the same procedure as before, and we obtain: 



-(1,1) 



(6.4) 



This is formally the same result as in the CDM case (see Appendix B for the leading terms of 7g\£ and T^g ), with 

the only difference that fs is replaced by fs- 

Let us now move to the evaluation of the leading contributions present in the genuine second-order part l2- The 
final results of Sect. V, for the particular case of a CDM model, can be easily generalized to the ACDM case 
starting from the observation that the leading contributions can only arise from terms containing the operators O^^ 



f2) z 

and Oj* in Eqs. (4.4) and (4.5). As a consequence, the first two terms 7^ and T2^' will give a subleading overall 

(2) (2) 

contribution, while the general result that the terms Tjj and give identically zero still holds 
leading contributions can be easily obtained, using the results in Eqs. (4.15) and (4.16), as follows: 



-(2) 



c(r55) 



1 



[BM - B3(77s)) + (54(770) - B4(?7.)) 



^3 

A r 



dii' (\B^{Tl') + Bi{i^')]k 



(6.5) 
(6.6) 
(6.7) 



These leading contributions can be now evaluated using Eqs. (4.8) and (4.9) and moving to redshift space, where we 
can write: 



Ho 



1 



1/2 



a„o(l + z)3 + l]AO ' 



where the suffix "0" appended to Vl„i and JIa denotes the present value of those fractions of critical density. We note 
that, as in the CDM case, the leading genuine second-order contributions are about one to two orders of magnitude 
smaller than the leading squared first-order contributions, evaluated above. 
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FIG. 2: A comparison of the primordial inflationary spectrum (long-dashed curve) with the spectrum of the CDM model 
neglecting baryons (thick solid curve) and of a ACDM model (thin solid curves), at various values of z. The dotted curves for 
the ACDM case describe the spectrum obtained by neglecting the baryon contribution (hence without taking into account the 
Silk-damping effect). 



We need now to insert the explicit form of the power spectrum. Considering the general solution of Eq. (4.2) we 
can also re-express the z-dependence of the power spectrum, in the linear regime, as follows: 

V^{k^z)=(^^ ^\T\k){^^^ , (6.9) 

where the previous CDM result (5.17), based on the transfer function T(fc) given in [10], is modified by the presence 
of the factor g{z)/gao originating from the time dependence of the gravitational perturbations. Another modification 
with respect to the CDM result, implicitly contained into the transfer function T{k), concerns the different numerical 
value of the equiUbrium scale A:eq (which now turns out to be lower because of the lower value of fimo)- The effects 
of such modifications are illustrated in Fig. 2 by comparing the ACDM spectrum of Eq. (6.9), at different values of 
z, with two z-independent spectra: the primordial spectrum of scalar perturbations (3/5)^A^, and the "transferred" 
spectrum of the CDM model, introduced in Sect. V (notice that, as expected, the ACDM and CDM spectra tend 
to coincide at large enough values of z). In the ACDM case the solid curves are obtained with a transfer function 
which takes into account the presence of baryonic matter, while the dashed curves correspond to the transfer function 
T{k) = To{k) of Sect.VC (without baryons). Here, however, we are always using a spectrum evaluated in the linear 
regime. 

As illustrated in Fig. 2, the effect of neglecting the baryonic fraction of ftm (and thus the associated Silk-damping 
effect) may lead to an overestimation up to 40% of the corresponding transfer function for scalar perturbations, in 
the range k > O.Ol/i Mpc""'^ (see [10]). In order to take into account this baryonic contribution we have to replace the 
value of q used in the previous section, i.e. q — /c/(13.41fceq), with the more accurate value given by [10]: 

, = (6.10) 

^ = ^-K"^ + TT(^)' ^'-''^ 

ar = 1 - 0.328 ln(43ir2™o/i^) +0.38 ln(22.3r2™o/i^) , (6.12) 

i KnO \ ' 'mO / 

44.51n(9.83/l]„ofe^) ^^ 
s = — Mpc. 6.13 
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FIG. 3: The fractional correction to the flux of Eq. (3.3) (thin curves) is plotted together with the fractional correction to 
the luminosity distance fd of Eq. (3.5) (thick curves), for a ACDM model with f^Ao — 0.73. We have used two different cutoff 
values: kuv = O.lMpc"^ (dashed curves) and kuv = IMpc ^ (solid curves). The left panel shows the results obtained with 
a linear spectrum without baryon contributions. The right panel illustrates the effects of including baryons (we have used, in 
particular, Qto = 0.046). 

Here flbo is the baryon density parameter, s is the sound horizon and T is the fc-dependent effective shape parameter. 

We have compared the above transfer function to the one which includes baryon acoustic oscillations (BAO) [10], 
and to a transfer function calculated numerically by using the so-called "code for anisotropics in the microwave 
background" (CAMB) [30]. We have checked, in particular, that the above simple form of transfer function is 
accurate to within a few percent compared to the one calculated numerically by CAMB, for all scales of interest. In 
addition, the effect of including BAO only produces oscillations of the spectrum around the above value. Since we 
are considering here integrals over a large range of k, the presence of BAO has a negligible effect on our final results, 
and will be neglected in the rest of the paper. 

We are now in the position of computing the fractional corrections to the averaged flux variable in a perturbed 
ACDM geometry. As discussed before, there are complicated average integrals which can be performed only by using 
some approximations. Once this is done, the remaining integration over k can be done numerically, exactly as in the 
case of the CDM model. 

In a ACDM context we may generally expect smaller corrections to the averaged flux, due to the fact that the 
perturbation spectrum V,p is suppressed by the presence of g{z). In addition (and as already stressed) the transfer 
function [10] turns out to be suppressed, at large k, because of a smaller value of the parameter fcgq (see Eq. (5.19)). 
These expectations are fully confirmed by an explicit numerical computation of |/$|, which we have performed with 
and without the inclusion of the baryon contributions into the transfer function T(k). 

The results of such a computation are illustrated in Fig. 3, and a comparison with Fig. 1 clearly shows that |/$| 
is smaller in the ACDM case than in the CDM case, and further (slightly) depressed when we take into account 
the presence of a small fraction of baryon matter (the curves presented in the right panel). In any case, the small 
values of |/$| at relatively large z, for a realistic ACDM scenario, lead us to conclude that the averaged flux is a 
particularly appropriate quantity for extracting from the observational data the "true" cosmological parameters. As 
we will discuss now, the situation is somewhat different for other functions of d^. 

B. Second-order corrections to other observables and dispersions 

Let us now consider other observables, beyond the flux, to see how the impact of the inhomogeneities may change. 
We will treat, in particular, the two important examples introduced in Sect. Ill, namely the luminosity distance cLl 
and the distance modulus /x. The fractional corrections to their averages (see Eqs. (3.4-3.6)) are qualitatively different 
from those of the averaged flux (represented by /$), because of the presence of extra contributions, unavoidable for 
any non-linear function of the flux and proportional to the square of the flrst order fluctuation ($i/<I>o)^. 

For a better understanding of such contributions let us start with the results obtained in [8] , concerning the second- 
order perturbative expansion of the luminosity distance, cLl = d^^^ + + and summarized in Sect. II B (see in 
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particular Eqs .(2.13) and (2.15)). Using those results we obtain: 

-Ii + (t.d.)« = 2 ( -EsJ- 
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(6.14) 



In order to determine the leading corrections we first notice that, by applying the procedure of Sect. V, and computing 
the averages for the CDM case, we obtain, to leading order, 
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Working in the context of a ACDM model, considering only these leading terms, and limiting ourselves to the linear 
regime, we find that the terms multiplying on the right hand side of the above equation can be calculated without 
approximations (as seen in the previous subsection). Also, we find that their contribution is controlled by the 
spectral factor k'^V^{k,r]o)- The first term, on the contrary, is due to to the so-called "lensing effect", dominates at 
large z (as already shown in [9] for a CDM model), and has leading spectral contributions of the type k^V^{k,r]o)- 
For such term, however, the integrals over time cannot be factorized with respect to the k integrals, and we must 
use the approximation already introduced in the previous subsection (namely, we have to replace in the integrands 
i^iv'jVo — v'l^^'^) with il;{ris,rio — rj' ,9°')). The full result for the spectral coefficient can then be finally written as 
follows: 



C(($i/ci>o)i) ^ + fc^ + 
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Because of the new term (due to lensing) affecting the average of fi and ~ but not of the flux $ - we may expect 
larger fractional corrections for these variables, as well as for other functions of $, at higher redshifts. This is indeed 
confirmed by the plots presented Fig. 3 reporting the results of an explicit numerical integration and comparing, in 
particular, the value of fd with the absolute value of /$. We obtain |/$| <^ fd, at large values of z where the lensing 
term dominates, both in the presence and in the absence of the baryon contribution to the total energy density. It 
should be stressed, however, that also the new fc^-enhanced contributions are free from IR and UV divergences, at 
least for the class of models we are considering. 

Let us now discuss to what extent the enhanced corrections due to the square of the first-order flux fluctuation 
can affect the determination of the dark-energy parameters, if quantities other than the flux are used to fit the 
observational data. To this purpose we may consider the much used (average of the) distance modulus given in Eq. 
(3.6), referring it, as usual, to a homogeneous Milne model with = 5logiQ[{2 + z)z/{2Hq)]. Considering Eqs. (3.6) 
and (3.9), where the averaged value of the distance modulus and its dispersion are given in function of /$ and of 

($i/$o)^^j we have investigated the magnitude of the effect in this case. In particular, in Fig. 4 we have compared 

the averaged value (/i) — with the corresponding expression for homogeneous ACDM models with different values 
of r^Ao- We have also illustrated the expected dispersion around the averaged result, represented by the dispersion 
previously reported in Eq. (3.9) (and already computed in [9] for the CDM case). 

We have found that the given inhomogeneities, on the average, may affect the determination of ^Iao obtained from 
the measure of the distance modulus, at large z, only at the third decimal figure (at least if the spectral contributions 
are computed in the linear regime). As we can see from Fig. 4, the curves for (/i) and for the corresponding 
unperturbed value ^^lrw ^-^^rith the same Oaq) practically coincide at large enough z. It should be stressed, also, 
that the dispersion on the distance modulus computed from Eq. (3.9) reaches, at large redshift, a value which is 
comparable with a change of about 2% in the dark energy parameter fl/^Q (cf. [1], considering however that baryonic 
effects have now been added in the transfer function). We shall see in the next section that this effect is enhanced by 
the use of non-linear power spectra. 



VII. THE ACDM MODEL: POWER SPECTRUM IN THE NON-LINEAR REGIME 

The linear spectra considered so far are sufficiently accurate only up to scales of order 0.1ft. Mpc^^. If we want to 
better study the effect of shorter-scale inhomogeneities on our light-cone averages we need to go beyond such linear 
approximation, taking into consideration the non-linear evolution of the gravitational perturbations. This can be done 
by using the so-called "HaloFit" models, which are known to reproduce quite accurately the results of cosmological 
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FIG. 4: The averaged distance modulus (fj.) — p*^ of Eq. (3.6) (thick solid curve), and its dispersion of Eq. (3.9) (shaded 
region) are computed for Qao = 0.73 and compared with the homogeneous value for the unperturbed ACDM models with, 
from bottom to top, Qao = 0.69, 0.71, 0.73, 0.75, 0.77 (dashed curves). We have used kuv = lMpc~^. The left panel shows 
the results obtained with a linear spectrum without baryon contributions. The right panel illustrates the effects of including 
baryons, with Qto ~ 0.046. 



A^-body simulations. In particular, the HaloFit model of [11] and its recent upgrade of [12] provide an accurate fitting 
formula for the power spectrum up to wavenumber k ~ 30/iMpc~^. 

In our previous paper [1] the analysis was limited to fc < 1 Mpc~^, but with the use of the non- linear power spectra 
we can now extend our analysis up to the maximum scale of the mentioned HaloFit models. The main difficulty with 
such an extension is that the time (i.e. z) dependence of the spectrum becomes more involved than in the linear case, 
since different scales no longer evolve independently. Hence, the need for introducing approximations in performing 
the integrals becomes even more essential in this case. 

Let us start by recalling a few details of the HaloFit model. The fractional density variance per unit Infc is 
represented by the variable A^(fc), defined by [11]: 



a' ^ 6{x)6ix) = I ^\S~,\' = I A^{k) d\nk , 



(7.1) 



which implies 
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where 6{x) = Sp{x)/p is the fractional density perturbation of the gravitational sources, and Sk is the associated 
Fourier component. On the other hand, the power spectrum of scalar perturbations, V^{k,z), is related to A^(fc,z) 
by the Poisson equation, holding at both the linear (L) and non-linear (NL) level [27]: 



(7.3) 



The linear part of the spectrum is used to introduce the normalization equation, defining the non-linearity length 
scale k~^{z), as follows: 



Al{k,z) exp(-(fc/fc. 



d\nk= 1. 



(7.4) 



It is then obvious that the scale fcg. is redshift-dependent. Since the non-linear power spectra obtained from the models 
[11] and [12] are using such a scale, they will be characterized by an implicit z-dependence which is more complicated 
than the one following from the usual growth factor of Eq. (4.3), and which cannot be factorized. 

Besides ka, the linear spectrum also determines two additional parameters, important for the construction of the 
HaloFit model: the effective spectral index and the parameter C, controlling the curvature of the spectral index 
at the scale fco-. They are defined by: 
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k [h Mpc"^] 

FIG. 5: The linear spectrum V^p (dotted curves) and the non-Unear spectrum for the HaloFit model of [11] (dashed curves) 
and of [12] (solid curves). In all three cases the spectrum is multiplied by (for graphical convenience) and is given for z = 
(thin curves) and z = 1.5 (thick curves). We have included baryons with 0.^0 = 0.046. 

Once the values of fc^, Hoff and C are determined, they can be inserted into a given HaloFit model [11, 12] to produce 
the non-linear power spectrum A^L(fc, z). One finally goes back to V^^{k, z) using again the Poisson equation (7.3). 

The non-linear spectra obtained in this way, and based on the previous linear spectrum with baryon contributions, 
are illustrated in Fig. 5 for the two HaloFit models [11] and [12]. The results are compared with the spectrum of the 
linear regime, for different values of the redshift (z = and z — 1.5). We can see that the spectra intersect each other 
in the non-linear regime (fc > 0.1/iMpc"^), as a result of the intricate redshift dependence. We can also observe that 
the two HaloFit models lead to the same result at low values of k (including baryons, but without BAO). 

A. Numerical results and comparison with the linear regime 

In order to evaluate our integrals using the non-linear power spectra we need to face the additional (already 
mentioned) problem due to the fact that the non-linear spectra - unlike the linear ones - cannot be factorized as a 
function of k times a function of z. In that case the full two-dimensional integration is highly non-trivial, and we have 
thus exploited a further approximation. In the presence of time-integrals of the mode functions (like those appearing, 
for instance, in the leading terms of (Ii^i)), we have parametrized the non-linear power spectrum in a factorized form 
as follows: 

v^^ik,z)^4rlr^^ik,z*), (7.6) 

and we have chosen z* = Zs/2 to try to minimize the error. 

Let us briefly comment on the validity of the approximation we have introduced. Our analysis being focussed on 
the redshift range 0.015 < z < 2, the above approximation is most inadequate only when we consider Zg ~ 2 (i.e. 
z* = I), and our mode functions are evaluated inside the integrals at the values of z most distant from z* (namely, 
z = 0.015 and z = 2). In that case we are lead to underestimate the spectrum by about a 40% factor for z = 0.015, 
and to overestimate it by about a 80% factor for z = 2. However, this only occurs in two narrow bands of k centered 
around the values k = I/iMpc~^ and k = 2/iMpc~^, while, outside these bands, our approximation is good. Since 
these errors are limited to only a part of the region of integration, both in z and in k, we can estimate an overall 
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accuracy at the 10% level, at least, for the results given by the adopted approximation. We have also checked the 
sensitivity of the numerical results to changes in z*, such as z* = Zs, and checked that our final results are only weakly 
dependent (typically at the 1% level) on the choice of z* . 

Once we have established the range of validity of the parametrization (7.6), we proceed in evaluating the z (or rf) 
integrals as we did in Sect. VI. The final results of this procedure are illustrated by the curves plotted in Figs. 6 
and 7, computed with the non-linear power spectrum following from the HaloFit model of [12] and including baryon 
contributions^. As illustrated for instance in Fig. 6, the fractional correction to dj^ turns out to be of order of a few 
parts in 10"'^ around z = 2, and smaller in the rest of the intermediate redshift range relevant for cosmic acceleration. 
On the other hand, in the same redshift range, the fractional correction to the flux is about two orders of magnitude 
smaller. 

Comparing with the results obtained with the linear spectrum (see Figs. 3, 4), we can see that taking into account 
the non-linearity distortions (and using higher cut-off values) enhances the backreaction effects on the considered 
functions of the luminosity distance, but not enough to reach a (currently) observable level. On the other hand, the 
dispersion, already large in the linear case, is further enhanced when non-linearities are included. In particular, the 
dispersion on ^ due to inhomogeneities is of order 10% around z = 2, which implies that the predictions of homogeneous 
models with f^Ao ranging from 0.68 to 0.78 lie inside one standard deviation with respect to the averaged predictions 
of a perturbed (inhomogeneous) model with JIao — 0.73. 
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FIG. 6: The fractional correction to the flux (/$, thin curves) and to the luminosity distance (/d, thick curves), for a perturbed 
ACDM model with JIao = 0.73. Unlike in Fig. 3, we have taken into account the non-linear contributions to the power spectrum 
given by the HaloFit model of [12] (including baryons), and we have used the following cutoff values: kuv = lO/iMpc"^ (dashed 
curves) and kuv = 30/iMpc~^ (solid curves). 



B. Comparing theory and observations via the intrinsic dispersion of the data 

Let us now consider in more detail our prediction for the dispersion cr^ induced by the presence of the inflationary 
perturbation background, and compare it with the intrinsic dispersion of the distance modulus that can be inferred 
from SNe la data. Our results for the dispersion are already implicitly contained in Fig. 7 but, for the sake of 
clarity, we have separately plotted our value of cr^ in Fig. 8, where the thick solid curve represents the value of cr^ 
obtained from Eq. (3.9), and plotted as a function of z. We can see from the figure that cr^ has a characteristic 



^ The two HaloFit models quoted above give similar results, and we have chosen to present here, for simplicity, only those obtained with 
one of them. 
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FIG. 7: The averaged distance modulus {^) — fi^' of Eq. (3.6) (thick solid curve), and its dispersion of Eq. (3.9) (shaded 
region), for a perturbed ACDM model with f^Ao ~ 0.73. Unlike Fig. 4, we have taken into account the non-linear contributions 
to the power spectrum given by the HaloFit model of [12] (including baryons), and used the cut-off kuv ~ 30ftMpc~^. The 
averaged results are compared with the homogeneous values of predicted by unperturbed ACDM models with (from bottom 
to top) Qao ~ 0.68, 0.69 0.71, 0.73, 0.75, 0.77, 0.78 (dashed curves). The right panel simply provides a zoom of the same 
curves, plotted in the smaller redshift range 0.5 < z < 2. 

z-dependence, with a minimal value of about 0.016 reached around z = 0.285. Also, the total dispersion of Eq. 
(3.9) nicely interpolates between the leading Doppler contribution obtained from the two last terms of Eq. (6.15) 
(represented by the dashed curve approaching zero at large z), and the leading lensing contribution obtained from 
the first term of Eq. (6.15) (represented by the dashed curve approaching zero at small z). 

The total variance a'^^ associated with the observational data, on the other hand, can be decomposed in general 
as follows (see e.g. [31-33]): 

(0' = (^r)' + K)' + «)' ■ (7-7) 

Here cr^* is the statistical uncertainty due, for instance, to the method adopted for fitting the light curve (e.g. the 
so-called SALT-II method [34]), but also to the uncertainty in the modeling of the supernova process. The term 
represents instead the uncertainty in redshift due to the peculiar velocity of the supernova as well as to the precision 
of spectroscopic measurements. Finally, cr™* is an unknown phenomenological quantity, needed to account for the 
remaining dispersion of the data with respect to the chosen homogeneous model. This part of the dispersion can 
be subsequently redefined whenever we are able to estimate some of the possible contributions it contains. The 
contribution we are mainly interested in here is the one originating from the lensing effect, which is dominant at large 
redshift. We can thus write, at large z: 

= + «"^)' , (7.8) 

where cr™* is the remaining source of intrinsic dispersion. 

Given the typical precision of current data [34, 35], a reasonable fit of the Hubble diagram does not seem to require 
a strong z-dependence of the parameter crj,"': for instance, a nearby sample gives a™* = 0.15 ± 0.02, to be compared 
with the value cr'"* = 0.12 ± 0.02 obtained for distant supernovae. On the other hand, as illustrated in Fig. 8, the 
results of our computations at z > 0.3 is very well captured by a linear behaviour which can be roughly fitted by 
a^°'^^{z) = 0.056z. We should also note that this contribution stays below 0.12 up to z ^ 2, which makes it perfectly 
compatible with observations so far performed. 

It is remarkable that the (above mentioned) simple linear fit of our curve for cr^(z) at z > 0.3 turns out to be 
very close to the experimental estimate reported in [36], namely crjf"*' = (0.05 ± 0.022)z. Also, such a fit is well 
compatible with the results of [37], namely ajf"*' = (0.055^g Q4^)z. These two observational estimates of the lensing 
dispersion, with the relative error bands, are illustrated by the shaded areas of Fig. 8. Our result is also in relatively 
good agreement with the simulations carried out in [38] which predicted an effect of 0.088z. Other authors ([39-41]), 
however, have found no indication of this z-dependence of cr^, and it is true that such a signal waits for a better 
observational evidence. 

It is likely that future improvements in the accuracy of SNe data will detect (or disprove) this effect at a higher 
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FIG. 8: The z-dependence of the total dispersion is illustrated by the thick solid curve, and it is separated into its "Doppler" 
part (dashed curve dominant at low z) and "lensing" part (dashed curve dominant at large z). The slope of the dispersion in 
the lensing-dominated regime is compared with the experimental estimates of Kronborg et al. [36] (dark shaded area), and of 
Jonsson et al. [37] (light shaded area). 

confidence level. In this respect, our results for a^{z) stand out as a challenging prediction^, which, we believe, could 
represent a further significative test of the concordance model. 

VIII. CONCLUSIONS 

Starting from the result of a previous paper [8] , where the luminosity-rcdshift relation has been computed to second 
order in the Poisson gauge, wc have proceeded here to the evaluation of the effects of a realistic stochastic background 
of perturbations on the determination of dark energy parameters. The basic tool we have used is the gauge-invariant 
light-cone averaging procedure proposed in [3], applied to different functions of the luminosity distance cLl averaged 
over a constant-redshift surface lying on our past light-cone. 

As already explained in [9], different functions of di differ by the sensitivity of their light-cone averages to fluctua- 
tions. Remarkably, a directly observable variable like the luminosity flux, $ ~ turns out to be the least sensitive 
to perturbations. Its averaged expression is also the simplest and, fortunately enough, other averages are readily 
computed once the one of the flux is given. Similarly, calculation of the dispersion is straightforward. 

As far as modeling inhomogeneities goes, we have used a concordance ACDM model with an arbitrary ^Im), upon 
which we have added a realistic spectrum of stochastic scalar perturbations up to second order (and we have also 
explained why, to this order, we have no contribution from vector and tensor perturbations). The perturbations have 
been taken as those originating from a quasi-scale-invariant primordial spectrum with a realistic transfer function 
(including baryons and the relative Silk-damping) after it undergoes a non-linear evolution according to the so-called 
HaloFit model for structure formation. This model appears to agree well with numerical A^-body simulations as well 
as with large-scale structure data. 

Our main conclusions (already succinctly presented in [1] for the case of a perturbation spectrum computed in the 
linear regime) are that the effect of perturbations on the averaged flux are extremely small, typically of order 10~^ at 
z ^ 0{1). Thus the average flux stands out as an extremely safe observable for determining dark-energy parameters 



We remark, incidentally, that our prediction for the Doppler-related dispersion at small z is also consistent with previous findings [42] 
on the so-called Poissonian peculiar- velocity contribution to cr^. We have checked that our result for the Doppler contribution (see Fig. 
8) is well fitted by an inverse power law: 0-^(2;) ~ 0.003232"^ in very good agreement, shape-wise, with the corresponding result in [42]. 
The fact that our prediction is about a factor 1.7 larger than the one of [42] is due, presumably, to the use of somewhat different power 
spectra (linear or non-linear, with or without baryons). 
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using the simplest FLRW geometry. Such observable is also practically insensitive (see Fig. 6) to the short-distance 
behaviour of the power spectrum. Other variables (like the commonly used distance modulus) receive corrections that 
are typically from two to three orders of magnitude larger (at large values of the redshift z), but still small-enough for 
allowing dark-energy measurements at the percent level without invoking theoretical corrections. On the other hand, 
they are more sensitive to the chosen value of the UV cutoff, and to the corrections arising from the power spectrum 
in the non-linear regime. 

We find, however, that the predicted intrinsic dispersion (or scatter) of the data due to just stochastic inhomo- 
geneities -and not to other well-known sources of dispersion- is considerably larger than their effect on averages. They 
imply that data should fluctuate in a band which, at large redshifts, covers the FLRW luminosity curves corresponding 
to a spread in J7ao of nearly 10% (see Fig. 7) . For limited statistics this irreducible dispersion will limit the precision 
with which dark-energy parameters can be extracted from the data. Particularly interesting, at large redshift, is the 
scatter due to lensing. Such an effect has been observed and a linear phenomenological fit to cr^ has been proposed 
[36, 37] with a slope dtr^/dz ~ 0.05 but with large errors. Our theoretical prediction is well described (for the consid- 
ered range of z) by the linear behaviour 17^(2:) ^ 0.056z which is not only consistent with the above phenomenological 
fits but also provides an interesting test of the concordance model if and when a more precise determination of crjj^"'' 
will become available. Also at small redshifts our (Doppler-induced) scatter, obeying an approximate inverse power 
law (T^t(z) ^ 0.00323z~^, looks compatible with observations and with previous theoretical estimates [42]. As a result 
of both effects we find the intrinsic dispersion cr^[z) to have a minimum of about 0.016 ai z 0.285. 

In any case, our conclusion is that, when averaging is applied to a physical observable within a well defined 
gauge invariant formalism, not even a small fraction of cosmological constant can be simulated by a stochastic 
(i.e. statistically homogeneous and isotropic) background of inhomogeneities. The situation is obviously different 
if one is willing to consider a deviation from the "almost scale-invariant" primordial spectrum, or a deterministic 
inhomogeneous and/or anisotropic cosmological model where we move along a very special geodesic, or one is ready 
to depart from the General Relativity framework. 

An interesting property of our averaging procedure is that, unlike other (more formal) definitions [25, 43], it leads 
to results which, for any realistic inhomogeneity power spectrum, are free from IR as well as UV divergences. The 
former property is very likely related to the gauge invariance of our procedure which, by definition, is unaffected by 
gauge artifacts due to super-horizon scales. Insensitivity to the UV regime is consistent with the intuition that very 
short-scale inhomogeneities should average out when considering large-scale physical observables (it may fail, instead, 
for artificially-defined spatial averages). Yet, some sensitivity to the actual UV cutoff remains in quantities that are 
controlled by a high (e.g. 3rd) moment of the power spectrum. Although, strictly speaking, we have only checked 
these nice properties up to second order, we strongly believe that they persist at higher orders as well. Actually, 
the use of our particularly suitable geodesic light-cone gauge [3] for performing light-cone averages may allow for a 
non-perturbative treatment of the backreaction problem. 
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Appendix A. Second-order vector and tensor perturbations 

We have already stressed in Sect. II B that vector and tensor perturbations automatically appear, at second order, 
sourced by the squared first-order perturbation terms. Hence, vector and tensor perturbations must be included 
in a consistent second-order computation of the luminosity distance, even if their contributions is negligible at first 
order (as expected, in particular, for a background of super-horizon perturbations generated by a phase of slow-roll 
inflation) . 
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Working in the Poisson gauge, and moving to spherical coordinates = (r, 0, 4>), we can rewrite the relevant part 
of the PG metric (2.11) as follows: 

ds%Q = [-drj^ + 2vid-qdx''] + [(7o)ij + Xij] dx'^dx^ , (A.l) 

so that 

^/ ) , (A.2) 

where 7^ = diag(l, r^^, r^^ sin^^ 9), and where we have called and x*"* the vector and tensor perturbations written 
in spherical polar coordinates. They satisfy the conditions ViU* = = V^x*^ and Xij = 0, where V; is the covariant 
gradient of three-dimensional Euclidean space in spherical coordinates. 

Following the same procedure as in the scalar case we can now evaluate the vector and tensor contributions to the 
coordinate transformation connecting Poisson and GLC gauge, and then express the perturbed GLC metric, up to 
second order, including the vector and tensor variables and x*"*- Such a detailed computation has already been 
performed, and its results presented in [8]. For the purpose of this paper it will be enough to recall here the vector 
and tensor contributions to the coordinate transformation between 6 and 9: 

r = r(o) + r(2) = r + ^' dx(^r{n+,x,e'^)-r''{v+.^.d'') + iS\v+.^.d'')J^dydha^T^+,y,e-)^ , (a.s) 

and to the 2x2 matrix 7°^ appearing in the GLC metric: 

/" dx9ji)''(r;+,a;,r)-x'''('7+,a;,r)+7o''(ry+,x,r) T 5da'-(r/+, y, r ) ) + (a O 6) 
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(A.4) 



where a'" = (w'^/2) — (x'"'^/4). Both results are needed, in fact, for the computation of the averaged flux (2.9). 

If we take into account v"^ and x*-* , and compute according to Eqs. (2.4) and (2.5), we find that the right-hand 
side of Eq. (2.13) has to be modified by the addition of a new term, (5y^(zs,0°), representing the effect of the vector 
and tensor part of the perturbed geometry (see [8] for its explicit expression). No modification is induced, however, 
on the corresponding equation for I(j,{zs) controlling the light-cone average of d^^, so that Eq. (2.14) holds even in 
the presence of vector and tensor perturbations. 

The sought average, in fact, is proportional to the proper area of the deformed two-sphere Ti{wo, Zs), and is given by 
I<l> ^ f ^Q\fl (see Eq. (2.9)). Considering the vector and tensor contributions to 7"-^ = det7'''' (obtained from Eq. 
(A.4)), and computing from Eq. (A. 3) the Jacobian determinant \dBldQ\^ we can express 1^ as an angular integral 
over the two-sphere with unperturbed measure d^il = siiiOdOdip. In that case many terms cancel among each other, 
and we end up with the result: 

I^{wo,Zs)-l = — I sin Ode I d<j>f{T],r,0,(P), (A.5) 
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where the integrand f{'q^r^9,<j)) is a simple expression proportional to the components of the vector and tensor 
perturbations. 

The above angular integrals are all identically vanishing, as we can check by expanding the perturbations in Fourier 
modes v\. and x^ ■ For each mode k we can choose, without loss of generality, the x'^ axis of our coordinate system 
aligned along the direction of k. Considering, for instance, tensor perturbations, we can then write the most general 
perturbed line-element, in Cartesian coordinates (omitting, for simplicity, the Fourier index), as follows: 

h^jdx'-dx^ = h+{ri,x^){dx^dx^ - dx^dx^) -f 2h-^{7],x'^)dx^dx'^ (A. 6) 

(we have called /i+ and h^, as usual, the two independent polarization modes). After transforming to spherical 
coordinates, using the standard definitions x^ — rsin^cosc/), x^ = rsinflsini/), x^ — rcosO, we easily obtain: 

sin ^20 1 

X'''' = sin^ 8{h+cos2(t) + hx sin 2(t>) ; x''" = (^+ cos 2(j) + hx sin 20) ; x''"^ = - (-^+ sin 2(t) + cos 20) ; 

2r r 

cos^ cos 

X^" = {h+ COS 24> + hy, sin 20) ; x""^ = ^ {-h+ sin 20 -f /ix cos 20) : 



r2 

1 



sin'^ ( 



(/i+cos20-H/ixSin20). (A. 7) 
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Since h+ = h+{r], r cos 9), hx = hx {1], r cos 9), all perturbation components depend on <j) only through cos 2(j) or sin 2(j>, 
so that their contribution averages to zero when inserted into Eq. (A. 5). The same is true for the case of vector 
perturbations, with the only difference that the </) dependence of v\ in spherical coordinates, is through cos (f> or sin0 
(corresponding to waves of helicity one instead of helicity two as in the tensor case). Also the vector contribution 
thus averages to zero when inserted into Eq. (A. 5). 



Appendix B. Computation of the C(7^'^^'^') spectral coefficients of 



We give here the result for the C(7^*-^'^'') spectral coefficients of (Xi^i) computed in the CDM case. We have 
introduced the convenient notation I = kArj, and we have enclosed in a box the leading contributions. Finally, we 
have defined Sinc(Z) = sin(/)/L 



Ar;2 3 



2S, 



Ar;2 3 



4S, 



+ {2 - 3^2 + (^2 _ 2) cos(/) + / sin(0 + ;3SinInt(/)} 



-4Ss4^(cos; - 2 + Sinc(Z) + /Sinlntf/), 
Arj 



'Arj^ 3 



CiT^'''^) = 4S,[1 - Sinc(0] - 4S,^ (cos(0 - Sinc(O) , 
3A(cos(0 - Sinc(Z)) - A_A_[2cos(/) + (l^ - 2)Sinc(/)] 



J fs ,2 



2S. 



fs ^^\^3 {\n^Avf] + n-iUoAr^ + - 6)] cos(Z) - 3(2 + HoA77)(-2 + ?')Sinc(0} , 

I. 



Us '^moAr^ 

, C{T^,'^)^Q,C{T^l''^) 





1 


■ 1 


1" 




" 3I2 + 


3I2 


" 6 



cos(0 - -Sinc(0 - -Sinlnt(/) 



^2 



1 



' UsA-n 

At]^ 3 



Vofs 

A772 



2 [1 - Sinc(0] + ^' [cos{l) - Sinc(0] 



At] 



[2cos(0 + (Z^-2)Sinc(0] 



4 <{ 1 - Sinc(/) + ^ [cos(0 - Sinc(0] \ , 



= 2S, 



Jo Is 

. Arj 

fo+fs 



(cos / - Sinc(O) + 



fo ;2 fofs 



3A?72 A7/2 



(2 cos(0 + {f - 2)Sinc(0) + ( 1 + 2 -^°^^^' ) (1 - Sinc(O) 



4 S.,^^^^ {1 - Sinc(0} - ySinInt(0 , C(r/6 ''^) - ySinInt(0 , C{r}r'^'') = ^ {-1 + cos(0 + I Sinlnt(0} 



UsAri 



2 - cos(/) - Sinc(/) - I Sinlnt(/) + ^ ( 1 



cos(Z) - I Sinlnt(0 



Vo 

At] 



[1 - Sinc(0] 



3/^ 



{2 - 3P + {f - 2) cos(0 + I sin(0 + /3SinInt(0} 



•H.Ary 



cos(0 + Sinc(/) + 



3A77 



A 
Ar] 



[2cos(0 + {l^ - 2)Sinc(0] } , CiT^T >) = 4 [1 - Sinc(0] 



4 + [-2(2 + 3f) + (4 + /2) cos(0 + I sin(0 + /(6 + P)SinInt(0] , C(r4' ^'') = 



3/2 
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